Investigating direct current potentials that affect native protein conformation during trapped ion mobility spectrometry-mass spectrometry.

Trapped ion mobility spectrometry-time-of-flight mass spectrometry (TIMS-TOFMS) has emerged as a tool to study protein conformational states. In TIMS, gas-phase ions are guided across the IM stages by applying direct current (DC) potentials (D1-6), which, however, might induce changes in protein structures through collisional activation. To define conditions for native protein analysis, we evaluated the influence of these DC potentials using the metalloenzyme bovine carbonic anhydrase (BCA) as primary test compound. The variation of DC potentials did not change BCA-ion charge and heme content but affected (relative) charge-state intensities and adduct retention. Constructed extracted-ion mobilograms and corresponding collisional cross-section (CCS) profiles gave useful insights in (alterations of) protein conformational state. For BCA, the D3 and D6 potential (which are applied between the deflection transfer and funnel 1 [F1] and the accumulation exit and the start of the ramp, respectively) had most profound effects, showing multimodal CCS distributions at higher potentials indicating gradual unfolding. The other DC potentials only marginally altered the CCS profiles of BCA. To allow for more general conclusions, five additional proteins of diverse molecular weight and conformational stability were analyzed, and for the main protein charge states, CCS profiles were constructed. Principal component analysis (PCA) of the obtained data showed that D1 and D3 exhibit the highest degree of correlation with the ratio of folded and unfolded protein (F/U) as extracted from the mobilograms obtained per set D potential. The correlation of D6 with F/U and protein charge were similar, and D2, D4, and D5 showed an inverse correlation with F/U but were correlated with protein charge. Although DC boundary values for induced conformational changes appeared protein dependent, a set of DC values could be determined, which assured native analysis of most proteins.

Structural analysis of ATP bound to the F1-ATPase ß-subunit monomer by solid-state NMR- insight into the hydrolysis mechanism in F1.

ATP-hydrolysis-associated conformational change of the ß-subunit during the rotation of F1-ATPase (F1) has been discussed using cryo-electron microscopy (cryo-EM). Since it is worthwhile to further investigate the conformation of ATP at the catalytic subunit through an alternative approach, the structure of ATP bound to the F1ß-subunit monomer (ß) was analyzed by solid-state NMR. The adenosine conformation of ATP-ß was similar to that of ATP analog in F1 crystal structures. 31P chemical shift analysis showed that the Pα and Pß conformations of ATP-ß are gauche-trans and trans-trans, respectively. The triphosphate chain is more extended in ATP-ß than in ATP analog in F1 crystals. This appears to be in the state just before ATP hydrolysis. Furthermore, the ATP-ß conformation is known to be more closed than the closed form in F1 crystal structures. In view of the cryo-EM results, ATP-ß would be a model of the most closed ß-subunit with ATP ready for hydrolysis in the hydrolysis stroke of the F1 rotation.

Terms and conditions.

Brian Plosky provides some context for a debate over the use of "intrinsically disordered" to describe regions of proteins.

Intrinsic disorder: A term to define the specific physicochemical characteristic of protein conformational heterogeneity.

In his commentary in this issue of Molecular Cell,1 Struhl reasons that the term "intrinsically disordered regions" represents a vague and confusing concept for protein function. However, the term "intrinsically disordered" highlights the important physicochemical characteristic of conformational heterogeneity. Thus, "intrinsically disordered" is the counterpart to the term "folded, " with neither term having specific functional implications.

Atomic force microscopy 3D structural reconstruction of individual particles in the study of amyloid protein assemblies.

Recent developments in atomic force microscopy (AFM) image analysis have made three-dimensional (3D) structural reconstruction of individual particles observed on 2D AFM height images a reality. Here, we review the emerging contact point reconstruction AFM (CPR-AFM) methodology and its application in 3D reconstruction of individual helical amyloid filaments in the context of the challenges presented by the structural analysis of highly polymorphous and heterogeneous amyloid protein structures. How individual particle-level structural analysis can contribute to resolving the amyloid polymorph structure-function relationships, the environmental triggers leading to protein misfolding and aggregation into amyloid species, the influences by the conditions or minor fluctuations in the initial monomeric protein structure on the speed of amyloid fibril formation, and the extent of the different types of amyloid species that can be formed, are discussed. Future perspectives in the capabilities of AFM-based 3D structural reconstruction methodology exploiting synergies with other recent AFM technology advances are also discussed to highlight the potential of AFM as an emergent general, accessible and multimodal structural biology tool for the analysis of individual biomolecules.

Structure-Based Protein Assembly Simulations Including Various Binding Sites and Conformations.

Many biological functions are mediated by large complexes formed by multiple proteins and other cellular macromolecules. Recent progress in experimental structure determination, as well as in integrative modeling and protein structure prediction using deep learning approaches, has resulted in a rapid increase in the number of solved multiprotein assemblies. However, the assembly process of large complexes from their components is much less well-studied. We introduce a rapid computational structure-based (SB) model, GoCa, that allows to follow the assembly process of large multiprotein complexes based on a known native structure. Beyond existing SB Go̅-type models, it distinguishes between intra- and intersubunit interactions, allowing us to include coupled folding and binding. It accounts automatically for the permutation of identical subunits in a complex and allows the definition of multiple minima (native) structures in the case of proteins that undergo global transitions during assembly. The model is successfully tested on several multiprotein complexes. The source code of the GoCa program including a tutorial is publicly available on Github: https://github.com/ZachariasLab/GoCa. We also provide a web source that allows users to quickly generate the necessary input files for a GoCa simulation: https://goca.t38webservices.nat.tum.de.

Functional protein dynamics in a crystal.

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

ALDELE: All-Purpose Deep Learning Toolkits for Predicting the Biocatalytic Activities of Enzymes.

Rapidly predicting enzyme properties for catalyzing specific substrates is essential for identifying potential enzymes for industrial transformations. The demand for sustainable production of valuable industry chemicals utilizing biological resources raised a pressing need to speed up biocatalyst screening using machine learning techniques. In this research, we developed an all-purpose deep-learning-based multiple-toolkit (ALDELE) workflow for screening enzyme catalysts. ALDELE incorporates both structural and sequence representations of proteins, alongside representations of ligands by subgraphs and overall physicochemical properties. Comprehensive evaluation demonstrated that ALDELE can predict the catalytic activities of enzymes, and particularly, it identifies residue-based hotspots to guide enzyme engineering and generates substrate heat maps to explore the substrate scope for a given biocatalyst. Moreover, our models notably match empirical data, reinforcing the practicality and reliability of our approach through the alignment with confirmed mutation sites. ALDELE offers a facile and comprehensive solution by integrating different toolkits tailored for different purposes at affordable computational cost and therefore would be valuable to speed up the discovery of new functional enzymes for their exploitation by the industry.

PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses.

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.

Main role of fractal-like nature of conformational space in subdiffusion in proteins.

Protein dynamics involves a myriad of mechanical movements happening at different time and space scales, which make it highly complex. One of the less understood features of protein dynamics is subdiffusivity, defined as sublinear dependence between displacement and time. Here, we use all-atoms molecular dynamics (MD) simulations to directly interrogate an already well-established theory and demonstrate that subdiffusivity arises from the fractal nature of the network of metastable conformations over which the dynamics, thought of as a diffusion process, takes place.

Intrinsically disordered regions (IDRs): A vague and confusing concept for protein function.

The term "intrinsically disordered region" (IDR) in proteins has been used in numerous publications. However, most proteins contain IDRs, the term refers to very different types of structures and functions, and many IDRs become structured upon interaction with other biomolecules. Thus, IDR is an unnecessary, vague, and ultimately confusing concept.

Conformational free-energy landscapes of a Na+/Ca2+ exchanger explain its alternating-access mechanism and functional specificity.

Secondary-active transporters catalyze the movement of myriad substances across all cellular membranes, typically against opposing concentration gradients, and without consuming any ATP. To do so, these proteins employ an intriguing structural mechanism evolved to be activated only upon recognition or release of the transported species. We examine this self-regulated mechanism using a homolog of the cardiac Na+/Ca2+ exchanger as a model system. Using advanced computer simulations, we map out the complete functional cycle of this transporter, including unknown conformations that we validate against existing experimental data. Calculated free-energy landscapes reveal why this transporter functions as an antiporter rather than a symporter, why it specifically exchanges Na+ and Ca2+, and why the stoichiometry of this exchange is exactly 3:1. We also rationalize why the protein does not exchange H+ for either Ca2+ or Na+, despite being able to bind H+ and its high similarity with H+/Ca2+ exchangers. Interestingly, the nature of this transporter is not explained by its primary structural states, known as inward- and outward-open conformations; instead, the defining factor is the feasibility of conformational intermediates between those states, wherein access pathways leading to the substrate binding sites become simultaneously occluded from both sides of the membrane. This analysis offers a physically coherent, broadly transferable route to understand the emergence of function from structure among secondary-active membrane transporters.

Multilevel superposition for deciphering the conformational variability of protein ensembles.

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

PURA syndrome-causing mutations impair PUR-domain integrity and affect P-body association.

Mutations in the human PURA gene cause the neurodevelopmental PURA syndrome. In contrast to several other monogenetic disorders, almost all reported mutations in this nucleic acid-binding protein result in the full disease penetrance. In this study, we observed that patient mutations across PURA impair its previously reported co-localization with processing bodies. These mutations either destroyed the folding integrity, RNA binding, or dimerization of PURA. We also solved the crystal structures of the N- and C-terminal PUR domains of human PURA and combined them with molecular dynamics simulations and nuclear magnetic resonance measurements. The observed unusually high dynamics and structural promiscuity of PURA indicated that this protein is particularly susceptible to mutations impairing its structural integrity. It offers an explanation why even conservative mutations across PURA result in the full penetrance of symptoms in patients with PURA syndrome.

PURA syndrome is a neurodevelopmental disorder that affects about 650 patients worldwide, resulting in a range of symptoms including neurodevelopmental delays, intellectual disability, muscle weakness, seizures, and eating difficulties. The condition is caused by a mutated gene that codes for a protein called PURA. PURA binds RNA ­ the molecule that carries genetic information so it can be translated into proteins ­ and has roles in regulating the production of new proteins. Contrary to other conditions that result from mutations in a single gene, PURA syndrome patients show 'high penetrance', meaning almost every reported mutation in the gene leads to symptoms. Proske, Janowski et al. wanted to understand the molecular basis for this high penetrance. To find out more, the researchers first examined how patient mutations affected the location of the PURA in the cell, using human cells grown in the laboratory. Normally, PURA travels to P-bodies, which are groupings of RNA and proteins involved in regulating which genes get translated into proteins. The researchers found that in cells carrying PURA syndrome mutations, PURA failed to move adequately to P-bodies. To find out how this 'mislocalization' might happen, Proske, Janowski et al. tested how different mutations affected the three-dimensional folding of PURA. These analyses showed that the mutations impair the protein's folding and thereby disrupt PURA's ability to bind RNA, which may explain why mutant PURA cannot localize correctly. Proske, Janowski et al. describe the molecular abnormalities of PURA underlying this disorder and show how molecular analysis of patient mutations can reveal the mechanisms of a disease at the cell level. The results show that the impact of mutations on the structural integrity of the protein, which affects its ability to bind RNA, are likely key to the symptoms of the syndrome. Additionally, their approach used establishes a way to predict and test mutations that will cause PURA syndrome. This may help to develop diagnostic tools for this condition.

Determining Protein Structures Using X-Ray Crystallography.

X-ray crystallography is a robust and widely used technique that facilitates the three-dimensional structure determination of proteins at an atomic scale. This methodology entails the growth of protein crystals under controlled conditions followed by their exposure to X-ray beams and the subsequent analysis of the resulting diffraction patterns via computational tools to determine the three-dimensional architecture of the protein. However, achieving high-resolution structures through X-ray crystallography can be quite challenging due to complexities associated with protein purity, crystallization efficiency, and crystal quality.In this chapter, we provide a detailed overview of the gene to structure determination pipeline used in X-ray crystallography, a crucial tool for understanding protein structures. The chapter covers the steps in protein crystallization, along with the processes of data collection, processing, structure determination, and refinement. The most commonly faced challenges throughout this procedure are also addressed. Finally, the importance of standardized protocols for reproducibility and accuracy is emphasized, as they are crucial for advancing the understanding of protein structure and function.

HYPK: A marginally disordered protein sensitive to charge decoration.

Intrinsically disordered proteins (IDPs) that lie close to the empirical boundary separating IDPs and folded proteins in Uversky's charge-hydropathy plot may behave as "marginal IDPs" and sensitively switch conformation upon changes in environment (temperature, crowding, and charge screening), sequence, or both. In our search for such a marginal IDP, we selected Huntingtin-interacting protein K (HYPK) near that boundary as a candidate; PKIα, also near that boundary, has lower secondary structure propensity; and Crk1, just across the boundary on the folded side, has higher secondary structure propensity. We used a qualitative Förster resonance energy transfer-based assay together with circular dichroism to simultaneously probe global and local conformation. HYPK shows several unique features indicating marginality: a cooperative transition in end-to-end distance with temperature, like Crk1 and folded proteins, but unlike PKIα; enhanced secondary structure upon crowding, in contrast to Crk1 and PKIα; and a cross-over from salt-induced expansion to compaction at high temperature, likely due to a structure-to-disorder transition not seen in Crk1 and PKIα. We then tested HYPK's sensitivity to charge patterning by designing charge-flipped variants including two specific sequences with identical amino acid composition that markedly differ in their predicted size and response to salt. The experimentally observed trends, also including mutants of PKIα, verify the predictions from sequence charge decoration metrics. Marginal proteins like HYPK show features of both folded and disordered proteins that make them sensitive to physicochemical perturbations and structural control by charge patterning.

Catalytic specificity and crystal structure of cystathionine γ-lyase from Pseudomonas aeruginosa.

The escalating drug resistance among microorganisms underscores the urgent need for innovative therapeutic strategies and a comprehensive understanding of bacteria's defense mechanisms against oxidative stress and antibiotics. Among the recently discovered barriers, the endogenous production of hydrogen sulfide (H2S) via the reverse transsulfuration pathway, emerges as a noteworthy factor. In this study, we have explored the catalytic capabilities and crystal structure of cystathionine γ-lyase from Pseudomonas aeruginosa (PaCGL), a multidrug-opportunistic pathogen chiefly responsible for nosocomial infections. In addition to a canonical L-cystathionine hydrolysis, PaCGL efficiently catalyzes the production of H2S using L-cysteine and/or L-homocysteine as alternative substrates. Comparative analysis with the human enzyme and counterparts from other pathogens revealed distinct structural features within the primary enzyme cavities. Specifically, a distinctly folded entrance loop could potentially modulate the access of substrates and/or inhibitors to the catalytic site. Our findings offer significant insights into the structural evolution of CGL enzymes across different pathogens and provide novel opportunities for developing specific inhibitors targeting PaCGL.

Can Protein Structure Prediction Methods Capture Alternative Conformations of Membrane Transporters?

Understanding the conformational dynamics of proteins, such as the inward-facing (IF) and outward-facing (OF) transition observed in transporters, is vital for elucidating their functional mechanisms. Despite significant advances in protein structure prediction (PSP) over the past three decades, most efforts have been focused on single-state prediction, leaving multistate or alternative conformation prediction (ACP) relatively unexplored. This discrepancy has led to the development of highly accurate PSP methods such as AlphaFold, yet their capabilities for ACP remain limited. To investigate the performance of current PSP methods in ACP, we curated a data set, named IOMemP, consisting of 32 experimentally determined high-resolution IF and OF structures of 16 membrane proteins with substantial conformational changes. We benchmarked 12 representative PSP methods, along with two recent multistate methods based on AlphaFold, against this data set. Our findings reveal a remarkably consistent preference for specific states across various PSP methods. We elucidated how coevolution information in MSAs influences state preference. Moreover, we showed that AlphaFold, when excluding coevolution information, estimated similar energies between the experimental IF and OF conformations, indicating that the energy model learned by AlphaFold is not biased toward any particular state. Our IOMemP data set and benchmark results are anticipated to advance the development of robust ACP methods.

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes.

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

Computational Insights into the Conformational Dynamics of HIV-1 Vpr in a Lipid Bilayer for Ion Channel Modeling.

HIV-1 Vpr is a multifunctional accessory protein consisting of 96 amino acids that play a critical role in viral pathogenesis. Among its diverse range of activities, Vpr can create a cation-selective ion channel within the plasma membrane. However, the oligomeric state of this channel has not yet been elucidated. In this study, we investigated the conformational dynamics of Vpr helices to model the ion channel topology. First, we employed a series of multiscale simulations to investigate the specific structure of monomeric Vpr in a membrane model. During the lipid bilayer self-assembly coarse grain simulation, the C-terminal helix (residues 56-77) effectively formed the transmembrane region, while the N-terminal helix exhibited an amphipathic nature by associating horizontally with a single leaflet. All-atom molecular dynamics (MD) simulations of full-length Vpr inside a phospholipid bilayer show that the C-terminal helix remains very stable inside the bilayer core in a vertical orientation. Subsequently, using the predicted C-terminal helix orientation and conformation, various oligomeric states (ranging from tetramer to heptamer) possibly forming the Vpr ion channel were built and further evaluated. Among these models, the pentameric form exhibited consistent stability in MD simulations and displayed a compatible conformation for a water-assisted ion transport mechanism. This study provides structural insights into the ion channel activity of the Vpr protein and the foundation for developing therapeutics against HIV-1 Vpr-related conditions.